Two-colour QCD at finite fundamental quark-number density and related 
theories.* 

S. J. Hands, a J. B. Kogut, b t S. E. Morrison, a and D. K. Sinclair. * 

a Department of Physics, University of Wales Swansea, Singleton Park, Swansea SA2 8PP, UK 
department of Physics, University of Illinois, 1110 West Green Street, Urbana, IL 61801-3080, USA 
C HEP Division, Argonne National Laboratory 9700 South Cass Avenue, Argonne, IL 60439, USA 

We are simulating SU (2) Yang- Mills theory with four flavours of dynamical quarks in the fundamental repre- 
sentation of SU(2) 'colour' at finite chemical potential, fi for quark number, as a model for QCD at finite baryon 
number density. In particular we observe that for /i large enough this theory undergoes a phase transition to a 
state with a diquark condensate which breaks quark-number symmetry. In this phase we examine the spectrum 
of light scalar and pseudoscalar bosons and see evidence for the Goldstone boson associated with this spontaneous 
symmetry breaking. This theory is closely related to QCD at finite chemical potential for isospin, a theory which 
we are now studying for SU(3) colour. 



1. Introduction 

QCD at finite baryon number density describes 
nuclear matter including neutron stars and heavy 
nuclei. Nuclear matter at finite temperature ex- 
isted in the early universe and may be observed 
in relativistic heavy-ion collisions at CERN and 
RHIC. 

QCD at finite chemical potential, [i for quark 
number has a complex fermion determinant. 
Hence current simulation algorithms fail for this 
theory. Working explicitly at finite baryon num- 
ber density trades this complex determinant for 
a sign problem. 

We therefore turn to studying models which 
have some of the properties of QCD at finite /i. 
One such property suggested for QCD at finite 
(i is diquark condensation. It is this which led 
us to study 2-colour (SU(2)) QCD with quarks 
in the fundamental representation of SU(2) co i our , 
which has a real non-negative determinant (and 
pfaffian), permitting simulations. 

2-colour QCD with fundamental quarks is con- 
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fining and has a sensible meson spectrum but 
an unphysical "baryon" spectrum. At m = 
and fi = the usual SU(Nf) x SU(N f ) x U(l) 
chiral symmetry is enlarged to SU(2Nf). The 
"meson" multiplets are enlarged to include 2- 
quark "baryons" . Spontaneous breakdown of chi- 
ral symmetry gives SU(2Nf) — > Sp(2Nf) rather 
than the usual SU R (N f ) x SU L (N f ) x U v (l) -> 
Uy(Nf). For large enough chemical potential we 
expect a spontaneous breakdown of quark num- 
ber with a diquark condensate. This condensate 
is a colour singlet and has associated Goldstone 
bosons |Q. For true QCD the condensate is, of 
necessity, coloured, and breaks the gauge symme- 
try in the Higgs manner - colour superconductiv- 
ity. For a chiral perturbation theory analysis of 
2-colour QCD at finite // see [pi. 

The 2 flavour version of this theory can be rein- 
terpreted as 2-colour QCD at finite chemical po- 
tential for isospin. However, it is easy to see that 
we can treat true (3-colour) QCD at finite chem- 
ical potential for isospin since it has a real, non- 
negative determinant. Since nuclear matter has 
finite (negative) isospin density as well as finite 
baryon number density, this theory describes a 
surface in the phase diagram for nuclear matter. 
For a discussion of this theory including a chiral 
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perturbation theory analysis see 

2. Lattice 2-colour QCD at finite quark- 
number chemical potential 

The standard staggered fermion transcription 
of this theory is: 

S f=J2 \ X[P(J*) + MX + ^[X T T2X + XT2X 7 



Pfaffian 
where 



Ar 2 A 
-A T Ar 2 



det(A*A + A 2 ) (2) 



A = Ip{p) + m 



(3) 



(1) 



where the chemical potential li is introduced by 
multiplying links in the +t direction by e M and 
those in the —t direction by e _A1 . The diquark 
source term (Majorana mass term) is added to 
allow us to observe spontaneous breakdown of 
quark-number on a finite lattice. 

Integrating out the fermion fields gives: 



Note that the pfaffian is strictly positive, so 
that we can use the hybrid molecular dynam- 
ics method to simulate this theory using "noisy" 
fcrmions to take the square root, giving Nj = 4. 

When m = A = [i = 0, the chiral symmetry of 
the above action is expanded to U(2), from the 
[7(1) x U(l)v of true lattice QCD. This breaks 
spontaneously - U{2) — > U (l)y - yielding 3 Gold- 
stone bosons. For m = A = 0,^^Owe expect 
spontaneous breakdown of quark number and 2 
Goldstone bosons - a scalar diquark and a pseu- 
doscalar diquark. For A = 0, m ^ 0, ^ we 
have no spontaneous symmetry breaking for small 
li. For li large enough (fj, > m^/2 ?) we expect 
spontaneous breakdown of quark number and one 
Goldstone boson - a scalar diquark. (See jl) for a 
full discussion of these symmetries, and for early 
simulations of the 8 flavour theory at A = 0.) 

We are simulating this Nf = 4 theory on an 
8 4 lattice, measuring the chiral and diquark con- 
densates, the fermion number density, the Wil- 
son/Polyakov line, etc.. We are repeating these 
simulations on 12 3 x 24 lattices, where, in addi- 
tion, we are measuring all local scalar and pseu- 
doscalar meson and diquark propagators (con- 
nected and disconnected). 



3. Results 

We have preliminary results for a relatively 
large quark mass, m = 0.1, and A = 0.01,0.02 
(and for small fi). We have recently started 
simulations at a small quark mass, m — 0.025 
and A = 0.0025, 0.005 where we should observe a 
more complex spectrum of Goldstone and pseudo- 
Goldstone bosons. Finally we will simulate at 
772 = where there should be 3 Goldstone bosons 
for (x = and 2 Goldstone bosons for /j, > 0. All 
our simulations are done at j3 = 1.5 w (3 c (N t = 
4). 

Figure [l] shows the diquark condensate (x Tr 2X) 
as a function of chemical potential potential. For 
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Figure 1. Diquark condensate as a function of li 
on an 8 lattice. The arrow indicates li « m^/2. 



small fx there is no diquark condensate and quark- 
number is a good symmetry. At /i — ii c ~ m 7r /2 
there is a phase transition, above which there is 
a diquark condensate, and quark-number is spon- 
taneously broken. This condensate increases to a 
maximum and then falls towards zero at large li. 
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Figure || shows the chiral condensate, (xx) as 
a function of fi. The chiral condensate is approx- 
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Figure 3. Quark- number density as a function of 
\i on an 8 4 lattice. 



Figure 2. Chiral condensate as a function of /i on 
an 8 4 lattice. 



imately constant for fi < [i c . For fi > fj, c it falls, 
approaching zero for large fj,. 

Figure |^ shows the quark-number density as a 
function of \x. The quark-number density jo is 
zero for fi < [i c and increases for /i > fi c ap- 
proaching the saturation value jo = 2 for large /i. 

Finally, in figure |i| we show the masses of the 
pion and the scalar diquark state which is orthog- 
onal to the diquark condensate, as functions of 
/i. This is from the 12 3 x 24 lattice where we 
have less complete 'data'. The mass of the scalar 
diquark falls roughly as the expected — 2/i 
as n is increased towards /i c . For /j, > /i c it is 
small enough and its A dependence is such that 
one can believe that it would become massless as 
A — > 0. Thus the scalar diquark does appear to be 
the Goldstone boson associated with spontaneous 
quark-number breaking. The pion mass appears 



to remain constant for fi < fj, c , dropping towards 
zero as \x is increased beyond fi c . This mass spec- 
trum and the behaviour of the order parameters 
is in good agreement with the predictions from 
chiral perturbation theory for \i < 0.6. 

4. QCD at finite isospin density 

For QCD at finite chemical potential for isospin 
(7 3 ), the staggered quark action is 

S f = J2 {xWi^) + m] X + Axir 2 e X } (4) 

sites 

where the explicit symmetry breaking term pro- 
portional to A is needed to observe spontaneous 
isospin breaking on finite lattices; Ti are isospin 
matrices acting on the isodoublet x> an d e = 
( — \Y +y+z+t . Integrating out the fermion fields 
yields: 

det(^U + A 2 ) (5) 
where 

A = lp(n) + m (6) 
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Figure 4. Pion (tt) and scalar diquark (gg) masses 
as functions of /z. The straight line is mass = 
TOjr — 2/i. The arrow is at fi = m^/2. 



Note .A is only a 1 x 1 matrix in isospin space. 
This describes 8 fermion flavours, so we use hy- 
brid molecular dynamics with "noisy" fermions 
to perform the simulations, to reduce this to 2 
flavours. Note this determinant has exactly the 
same form as that for 2-colour QCD at finite 
quark number density, except now we can use 
SU(2>) colour. To get a sensible flavour inter- 
pretation for this 2-flavour theory in the contin- 
uum limit requires making the field redefinition 
X2 — * £5X2, where £5 is the analogue of 75 in 
-St/ (4) flavour space. 

5. Conclusions 

2-colour QCD at finite chemical potential for 
quark number undergoes a phase transition at 
M = Mc < m % /2. The high /i phase is charac- 
terized by spontaneous breakdown of quark num- 
ber precipitated by a diquark condensate, with a 
single Goldstone boson. Is there a second transi- 
tion at high jjL to a free-field phase? What can we 



learn from this model about diquark condensates 
in true (S*t/(3)) QCD? 

We need to simulate at smaller quark masses 
(including zero) where the competition between 
chiral and quark-number symmetry breaking 
should lead to a more complex spectrum of Gold- 
stone and pseudo-Goldstone bosons. A more ex- 
tensive analysis of the spectrum of scalar and 
pseudoscalar mesons and diquarks is called for. 

QCD at finite chemical potential, /i/, for 
isospin maps a surface in the phase diagram of 
nuclear matter which, in addition to having a fi- 
nite baryon number density, also has a finite (neg- 
ative) isospin density. A reinterpretation of 2- 
colour QCD at finite quark number density yields 
2-colour QCD at finite isospin density, and gives 
us a guide as to what to expect in the 3-colour 
case. At finite [ij, true (3-colour) QCD has a real 
positive determinant which allows us to simulate 
it. For large enough chemical potential we should 
get spontaneous breaking of the remaining t/(l) 
isospin symmetry generated by I3 with a charged 
pion condensate and one true Goldstone pion. 

Recently, the work of has been extended to 
finite temperature as well as fx by the Hiroshima 
group 0. In addition the eigenvalues of the Dirac 
matrix in jlj] have been studied by ||. Finally 
we would like to mention related work on gauge 
theories with adjoint fermions Q. 
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